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Abstract 

We demonstrate that a single reconnection of two quantum vortices can lead to creation of 
a cascade of vortex rings. Our analysis, motivated by the analytical solution in LIA, involves 
high-resolution Biot-Savart and Gross-Pitaevskii simulations. The latter showed that the rings 
cascade starts on the atomic scale, with rings diameters orders of magnitude smaller than the 
characteristic line spacing in the tangle. So created vortex rings may penetrate the tangle and 
annihilate on the boundaries. This provides an efficient mechanism of the vortex tangle decay in 
very low temperatures. 
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We consider quantum vortex lines of constant circulation k (for the superffuid 4 He k = 
h/rriHe — 9.97 x 10~ 4 cm 2 /s). The curve traced out by a vortex filament is specified in 
the parametric form s(£, t), with t and £ denoting respectively time and arc length. The 
vortex local velocity s(£, £), given by the Biot-Savart (BS) integral, in some cases can be 
approximated in terms of the, so called, localized induction approximation (LIA) retaining 
only the effects of the local vortex curvature pQ; which in nondimensional units reads 

s((,t)=s'xs" + as", (1) 

where the overdot and the prime denote the derivatives with respect to t and £ respectively, 
and a is the non-dimensional friction parameter. For a > 0, as showed by Lipniacki [2H3], Eq. 
[T] has four-dimensional class of self-similar solutions. When the initial vortex configuration 
consists of two half-lines with a common origin, the line motion is equivalent to a homothety 
transformation s(£, t) = s(l)\/i, with I := £/yi. In terms of curvature c(£, t) = K(l) and 
torsion r(£, t) = T(l) the self similar solution can be given in the implicit form 



1{K) = W— r 1 (2) 
V Ik ks/ln(K /k) + a{iq - k 2 ) 



a 



K' 

where K = K{1 = 0). In the limit of a — > this solution converges to the self-similar 
solution found by Buttke [I], and analyzed by Svistunov [5 J 



We reconstructed vortex lines given by the solution (|2j|3j) using the Frenet-Seret equations 
to show that for sufficiently small a (a < 0.44) and sufficiently small angle between the 
reconnecting lines 7 < 7i/>i(a) the resulting vortex line has two or more self-crossings, see 
Fig. [T]A.. Such solutions are not physical, but their existence suggests that the reconnection 
of two straight vortex lines at a sufficiently small angle may lead to a series of vortex self- 
reconnections and the creation of a cascade of vortex rings of growing diameter, as predicted 
at T = by Svistunov [5]. We confirmed creation of vortex rings cascades by performing 
high-resolution BS numerical simulations (following the numerical method proposed by Aarts 
[6]), starting from the configuration which arises shortly after the reconnection of two straight 
lines Fig. [T]B. To avoid singularity, the sharp corner in the initial configuration was replaced 




FIG. 1: A) Self-similar solution in the LI A approximation of one of the two vortices created by 
the reconnection of two straight vortex lines at the angle 7 = 5° for a = 0.01. The presented 
'unphysical' solution has two self-crossings. B) Generation of a vortex ring in a Biot-Savart law 
simulation due to the self-reconnection of the vortex filament created by the reconnection of two 
straight vortex lines for the same 7 and a as in A). The initial configuration is shown as a grey 
line. 
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FIG. 2: The critical angle ^lia{oi) below which the analytical solutions have at least two self 
crossings and the computed critical angle 755(a) below which vortex rings are generated in the 
Biot-Savart simulations. For a > 0.3 nonlocal interactions result in the collapse of a whole vortex 
rather than in the creation of vortex rings. 

by an arc of radius three time larger than the radius of the vortex core, a s m 1.3A. In 
the example simulation showed in Fig. [T]B, performed for a = 0.1 and 7 = 5°, we observe 
the creation of a vortex ring. In Fig. [2] we compare the critical angle ^lia{ol) below which 
analytical solutions have at least two self crossings and the estimated critical angle 755(0;) 
below which the vortex rings are generated in the BS simulations. Non-local interactions 
present in the BS simulations enable ring generation for even larger 7 than expected from the 
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LIA solutions (755(0) ~ 10.4° versus ^lia(Q) ~ 8.5°). In addition the non-local interactions 
cause two straight nearly antiparallel vortex lines to approach closely which allow them to 
undergo reconnection and initiate the vortex loops cascade, Fig. [3]A 

Since the BS dynamics may describe the vortex motion before and after the reconnection 
but not the reconnection event itself, we repeated the simulation of the reconnection of two 
nearly antiparallel vortices based on the Gross-Pitaevskii (GP) equation, implementing the 
Dufort-Frankel scheme described in [7]. As shown in Fig. [3j3 the first ring arising in the 
cascade has the atomic scale radius of the order of 3a s . This explains why the phenomenon 
was overlooked in the large scale BS simulations of the tangle in which the cusps arising 
after reconnections are replaced by the arcs of radii comparable to the characteristic radius 
of curvature of lines in the tangle. Interestingly, in one of the first reconnection studies 
based on the GP equation, Koplik and Levine [8] showed a transient configuration, possibly 
preceding the separation of a tiny vortex loop, which either quickly collapsed or left the 
small simulation box. 

In further analysis we focus on the zero-temperature limit (a = 0) and restrict ourselves 
to numerically more efficient BS simulations. It follows from the similarity of the subsequent 
vortex rings, and is confirmed by the numerical simulations (see Fig. |4|3), that their lengths 
li form a geometric sequence, i.e. 

^ = *(7)- (5) 
The simulations shown in Fig. indicate that 37(7) is a monotonically growing function 
of the reconnection angle 7 diverging as 7 — > jbs and may be well approximated as (ft (7) = 
1bs/(ibs ~ 7), see Fig. |4p. Correspondingly, the times of the subsequent ring separations 
also form a geometric sequence 

-7^ = Qtin), where ^(7) = qf (7). (6) 

The last equality, implies that the line-length l(t) which evaporates in the form of vortex 
rings during time t after the reconnection of two straight vortex lines is (in dimensional 
units) l(t) = 0(7) y/Tti, with numerically estimated 0(7) = 2.5 ± 0.3 for < 7 < 7^5- Having 
numerically determined (#(7) one can calculate the line- length which evaporates in the form 
of vortex rings "%2 i k ~ qil n /{qi — 1), where l n is the length of the last ring of the cascade. 
For an idealized reconnection of two straight vortex filaments in infinite volume the cascade 
would be infinite. However, in the realistic situation of the reconnections of curved vortex 

4 




FIG. 3: Biot-Savart A) and Gross-Pitaevskii B) simulations of the reconnection of two initially 
straight vortex filaments in pure superfluid (a = 0), inclined at the angle 7 = 5°. In the BS 
simulation the initial separation is of 2 x 10 3 a s ; the sharp corner arising at the reconnection is 
smoothed by the arc of radius equal to 20a s . In the GP simulation the isosurfaces of |<1?| 2 = 0.3 are 
shown. The initial line separation is equal Aa s ; the two arising rings have radii of the order of 3a s . 

lines in the tangle, the cascade of rings will be terminated when the angle between the two 
lines becomes larger than the critical angle 7^5. In Fig [5] we show the most unfavorable 
configuration in which the angle between the reconnecting lines grows fastest. Assuming 
that vortex lines with radius of curvature r reconnect at angle 7^5/2, we may give the 
lower bound (based on the most unfavorable configuration) of the evaporated line length 
as Uot = 27rr7£S'/360 o . The characteristic radius of curvature for the line in the tangle 
is r = c^L -1 / 2 , where c\ is the nondimensional coefficient of curvature, introduced and 
estimated by Schwarz [9]. 

To estimate the effective vortex tangle decay rate in the T — y limit we distinguish two 
regimes of 1) 'optically' transparent and of 2) 'optically' opaque tangle. For the first regime 
we assume that all created rings are annihilated at the boundaries, and as a result the average 
line length loss due to each reconnection is l to t = ^IbsCi 1 L^ 1 / 2 /360° ~ 0.09cq ^^Z 2 . For 
the optically opaque tangle only the smaller rings may reach the boundary and annihilate. 
The mean free path S of a ring of diameter d in the tangle of density L can be estimated 
as S = l/(dL). Assuming, for simplicity, that only rings having the mean free path longer 
than half of the vessel diameter D/2 (i.e. rings having diameter d < 2/(DL)) will reach 
the boundary, we calculate the maximum possible line-length loss associated with a single 
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FIG. 4: Rings generation after the reconnection of two straight vortex lines (BS simulations at 
T = 0). A) vortex configuration prior to the second ring separation for reconnection angles 2°, 5°, 
7°. B) The lengths li of the subsequent rings for the reconnection angles 7 = 1°, 2°, 3°, 5°, 7° are 
shown in logarithmic scale. C) Inverse ratio of lengths of the subsequent rings (l/<ft, cf. Eq. [5]) as 
a function 7. The critical reconnection angle "fss 1S marked. 




FIG. 5: The BS simulations at T = showing rings generation after the reconnection of two 
co-planar rings with radii of 10 3 a s at angle 7 = 5°. The rings cascade progresses until the angle 
between the two arms exceeds the critical value 7ssi three out of four created rings are shown. 

reconnection by taking l n = 2n/(DL) 

Now, if l* tot {^Bs) > hot (which implies ^bsc^L^^D < 360°) the tangle is optically transpar- 
ent. Otherwise the tangle is optically opaque. The reconnection frequency per unit volume 
is f r = ci(3L 5 ^ 2 [10]. Assuming the random distribution of the reconnection angles, the 
fraction of reconnections at which vortex cascades are created is fi — 7bs/180 o and thus we 
obtain the decay rate of an optically transparent tangle as 

^ = -Vfrhot ^ -0.01/3L 2 , (8) 
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and for an optically opaque tangle ("fBs c i l L 1 l 2 D < 360°) as 



dL „ 1 



dt 7b 5 



7c riBS 



Pfr / kotfa + / Ul) d l 



</ 7c 



_2^n , J 2M L^D\Y (9) 
D V V 360 ° c i J J 
where 7 C = 737772- As found already by Schwarz [TJ, even if the lines are initially inclined 
at larger angle they reorientate so that at the moment of reconnection they become almost 
antiparallel. As a result one may expect that p, is substantially larger than 7^5/180° and 
thus the decay rate in Eq. [8] should be regarded as a lower bound. The vortex tangle decay 
rate for the optically transparent tangle regime has the same form as the decay term in the 
classical Vinen equation, which in Schwarz's notation [9] reads [dL/dt]d ec = —ctc\fiL 2 , where 
c 2 is the coefficient of averaged squared curvature. According to Schwarz's simulations j9] 
for a steady state tangle ac\ = 0.12, for the smallest studied a = 0.01 (corresponding to 
T = 1.07 K). Since at low temperatures, a ~ T 5 we may expect that dissipation due to 
rings generation exceeds the friction force dissipation at temperature below T ~ 0.65K. 
Buttke's solution [1], implies that the wave number k of the Kelvin helical waves at £ = sjld 
decreases as 1/ V«, while the amplitude A increases as \/k£ . Since the energy per unit 
length radiated by sound (for dipole radiation, neglecting logarithmic terms, see (TTJ) is 
II = pn 5 k 6 A 2 /(64:7rc1) (where c s is sound velocity) the rate of energy dissipation and line 
length reduction decreases as t~ 2 . Consequently the line length is lost due to sound emission 
immediately each after reconnection event. The associated decay rate is dL/dt = —f r l s = 
—ci(3l s L 5 ^ 2 , where l s is line loss associates with single reconnection. Assuming that l s is of 
order of k/c s (which defines the spatial scale of the effective sound emission) we obtain that 
the Kelvin wave dissipation mechanism dominates over the dissipation via rings cascades 
generation when 

K Tl/2 

pci > 0.01, (10) 

where p is a constant. The analogous estimation was derived by Vinen [TTJ who found that 
Kelvin wave dissipation dominates over the mutual friction dissipation (which has the same 
form as the dissipation rings cascades) for 

> a, (11) 

c s 

which for L = 10 6 /cm 2 gives the critical temperature of T ~ 0.65K. Consequently, inequal- 
ities 10 11 imply that there are three regimes: i) high-temperature regime (a > kL 1 ^ 2 /c s ) 



when mutual friction dominates; ii) the low-temperature, dense tangle regime (a < kL 1 / 2 /c s , 
L > 0.01c s /pciK ~ 10 6 ) when Kelvin wave dissipation dominates and iii) the low- 
temperature, sparse tangle regime (a < kL 1//2 /c s , L < 10 6 ) when the most efficient mecha- 
nism of dissipation is the generation of vortex rings. 

The question whether T — > is a singular or regular limit of quantum turbulence is still 
not resolved. Experiments on quantum turbulence imply that some dissipative mechanisms 
persist even at the lowest attainable temperatures p2] • The mechanisms invoked to explain 
this dissipation in pure superfluid are all associated with vortex reconnection. Reconnec- 
tions cause the direct line loss [13] and they trigger Kelvin waves propagating along the 
reconnected vortices [5j [HI [15] . The nonlinear cascade of such waves may possibly transfer 
energy to the smallest scales where it is dissipated by the emission of phonons [13 [IB] . Fi- 
nally, as demonstrated here, on the basis of the fine-scale Biot-Savart and Gross-Pitaevskii 
simulations, vortex reconnections at sufficiently small angle (7 < 10.4°) lead to the creation 
of cascades of vortex rings of diameters starting from the atomic scale. Creation of rings 
cascades introduces the qualitative difference in the line length dissipation. Although, the 
distribution of wave numbers corresponding to subsequent rings and that of Kelvin waves 
arising at larger reconnection angles are similar, the time evolution of Kelvin waves and 
of rings are different; the wave vectors of Kelvin waves quickly decrease in time and thus 
the sound emission stops, while the curvatures of vortex rings grow in the course of energy 
dissipation thus increasing the dissipation rate. 

In conclusion, creation of vortex rings cascade provides a novel, efficient mechanism of 
line-length loss at very low temperatures which is not based on acoustic dissipation. Ring 
creation and emission is an old idea. However, starting from the hypothetical Feynman's 
scenario in which vortex rings were assumed to split into smaller rings it was always assumed 
that energy cascades from small to large wavelengths. The evaporation via vortex ring gen- 
eration has been considered by Barenghi and Samuels [17]. However, also in this mechanism 
the radii of the generated vortex rings are of the order of the characteristic vortex line cur- 
vature in the tangle. Such vortex rings may escape only from a small, localized packet of 
vorticity. Here we showed that the classical large-to-small-scale cascade can be bypassed 
when the reconnecting vortex lines are at least locally nearly antiparallel. A sequence of 
vortex rings is then produced with the diameters starting at the scale of the core size and 
increasing up to the scale of characteristic radius of curvature in the tangle. Opposite to 
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the normally envisaged cascade (sometimes called Richardson cascade) this is a mechanism 
of direct line-length transfer to the smallest spatial scales. 
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